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摘 要 : 全 天 相机 云图 在 天 文 观测 领域 具有 重要 价值 ,但 云图 所 含 噪声 会 造成 不 良 影响 。 


本 文 先后 分 析 了 星光 噪声 、 月 光 噪 声 、 日 光 噪声 、 


雷电 噪声 、 反 光 噪声 和 特殊 噪声 等 ， 并 总 


结 出 不 同类 型 噪声 的 分 布 、 区 域 大 小 等 方面 的 特点 ,最 后 结合 与 云层 的 相对 位 置 关系 提出 了 


边界 阔 值 法 、 时 间 轴 法 、 反 比例 线性 变换 等 8 种 图 像 处 理 去 噪 算法 ， 并 形成 一 个 完整 的 去 


噪 系统 。 通 过 实验 对 不 同 的 算法 加 以 比较 


, 得 出 不 同 噪声 情况 下 的 最 佳 去 噪 方案 。 


实验 结果 


表明 , 该 系统 可 以 有 效 去 除 常见 噪声 ， 基 本 保持 云层 和 天 空 的 信息 完整 性 ,可 简单 便捷 实现 


BAAR. 
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0 引 言 
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总 体 而 言 , 云图 


图 分 析 造 成 干扰 ,使 云图 平均 灰 度 值 变化 或 将 噪声 部 
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全 天 相机 云图 是 通过 全 天 空 成 像 技术 ， 通 过 全 
理 和 显示 24 小 时 的 天 空 状态 。 全 天 相机 拍摄 的 云图 是 角度 自 下 而 - 


图 去 噪 , 包括 


如 完备 DCT ^ 
广义 交叉 验证 (GCV) 


图 像 边 缘 清 


引入 噪声 去 品 算 法 的 研究 较为 成 熟 。 故 本 文 主要 研究 云图 
法 ， 后 文 所 提 噪 声 均 为 云图 固有 噪声 。 
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自动 、 全 色彩 天 空 成 像 系统 


， 提供 实时 处 


上 的 、 以 天 衬 为 背景 的 云 


图 "'"， 换 言 之 ， 全 天 相机 云图 是 一 种 地 基 云 图 。 这 也 就 导致 全 天 相机 云图 的 固有 噪声 主要 
来 自 于 天 空中 的 因素 。 目 前 ， 常 采用 RGB 三 通道 中 的 方法 进行 云 检测 及 太阳 去 除 ， 但 由 于 
空 背景 的 不 均匀 性 , RGB 三 通道 法 存在 一 定 误差 。R. Chauvin "等 人 通过 生成 晴空 图 的 方 


法 对 RGB 三 通道 法 进行 改进 ， 
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除去 太阳 噪声 。 但 该 方法 需 结合 云图 每 个 像素 的 像素 / 
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KAA PZA) 和 像素 /方位 角 (PAA)， 对 相机 的 校准 要 求 较 高 ， 成 本 高 且 易 产生 误差 。R. 
Román" ^^& A diti Y CPC (Camera Plus Ceilometer) 算法 ,可 检测 太阳 是 否 被 云层 遮挡 ， 提 
高 了 云 检 测 准确 度 ， 但 该 方法 没有 除去 太阳 噪声 。J.Alonsot| 等 人 提出 了 结合 辐射 数据 确定 
云图 太阳 位 置 的 方法 ， 但 该 方法 需要 结合 辐射 数据 ， 不 能 直接 处 理 云图 。 胡 树 贞 9 等 人 提 
T zu pH A, 但 该 方法 只 是 削弱 了 太阳 噪声 的 影响 ， 并 没有 完全 除去 太阳 噪声 。 魏 
诗 雅 [7 等 人 提出 了 时 间 分 割 法 和 差分 法 除去 月 光 噪 声 ， 但 这 两 种 方法 仅 能 处 理 灰 度 值 图 像 ， 
无 法 保留 云图 的 RGB 信息 。 总 体 而 言 ， 国 内 外 对 云图 固有 噪声 去 噪 算法 的 研究 多 为 去 除 或 
削弱 太 阳 和 月 亮 噪 声 ， 而 云图 中 还 存在 如 雷电 、 冰 蝇 等 干扰 项 。 并 且 ， 许多 去 品 算 法 需要 结 
合 辐射 中 、 温 度 M9 等 数据 ， 较 为 繁琐 ， 或 对 相机 等 设备 有 较 高 要 求 (1 ， 成 本 较 高 。 本 文通 
过 图 像 处 理 技术 ， 可 实现 多 类 型 噪声 去 除 工作 ， 简 单 便捷 ， 成 本 较 低 。 

本 文 的 主要 创新 点 有 三 : 第 一 是 在 经 过 大 量 的 人 工 阅 图 后 将 噪声 进行 合理 分 类 , 把 云图 
回 有 噪声 甄别 为 $ 种 常见 噪声 和 4 种 特殊 噪声 ,并 对 应 于 不 同 的 天 空 状态 ; 第 二 是 根据 不 同 
噪声 的 特点 进行 有 针对 性 处 理 ， 总 结 出 不 同 噪声 的 分 布 、 区 域 大 小 等 特点 ， 然 后 结合 噪声 与 
云层 相对 位 置 等 因素 有 和 针对 性 提出 边界 阔 值 法 、 时 间 轴 法 、 反 比例 线性 变换 等 8 种 去 噪 算 法 ， 
从 而 基本 实现 各 种 天 空 状态 下 的 云图 去 噪 工作。 第 三 是 本 文 提出 了 一 种 交互 式 多 类 型 噪声 去 
除 系 统 ， 可 减少 算法 执行 时 间 ， 并 选取 200 张 典型 云图 噪声 数据 进行 实验 分 析 。 实 验 结果 表 
W, 该 系统 可 以 有 效 去 除 常见 噪声 ， 基 本 保持 云层 和 天 空 信息 完整 性 ,简单 便捷 实现 云图 去 
WR 

本 文 的 数据 来 源 是 中 国 科学 院 国 家 天 文 台 提供 的 LAMOST 的 云图 数据 ， 共 计 有 2201 
张大 小 为 2592*1728 像素 的 彩色 云图 。 


1 云图 噪声 分 类 及 去 噪 算法 


本 文 将 云图 噪声 分 为 星光 噪声 、 月 光 噪 声 、 日 光 噪 声 、 雷 电 噪声 、 反 光 噪 声 和 特殊 噪声 。 
1.1 星光 噪声 
星光 噪声 的 灰 度 级 跨度 较 大 , 像素 点 小 ， 灰 度 值 比 天空 高 。 因 此 可 以 认定 星光 是 一 种 随 
机 出 现 的 白 点 ， 是 夜空 这 一 瞳 的 区 域 中 的 白色 像素 。 这 一 特征 与 椒盐 噪声 的 特征 相同 。 故 可 
用 中 值 滤波 去 品 算 法 。 此 外 , 根据 全 天 相机 连续 拍摄 的 特性 ， 可 以 采用 的 第 二 种 去 噪 算法 为 
时 间 轴 法 。 在 针对 该 类 噪声 的 处 理 中 ， 所 有 的 彩色 云图 通过 式 (1-1) 全 部 转 为 256 级 灰 度 图 
后 再 行 处 理 。 


Gray = 0.299 *R+0.587*G+0.114*B (1-1) 
式 中 Gray 为 像素 点 的 灰 度 值 ， 范 围 为 0-255 。 
1. 1.1 中 值 滤波 法 : 中 值 滤波 法 是 去 除 单 点 无 规律 分 布 噪声 的 有 效 手 段 。 该 方法 通过 确定 一 
个 含 奇 数 像素 的 滑动 窗口 , 将 窗口 内 的 各 像素 灰 度 值 排序 , 用 其 中 值 代 替 原 灰 度 值 成 为 窗口 
中 心 的 灰 度 值 ， 从 而 消除 孤立 的 星光 噪声 点 。 它 属于 一 种 非 线性 的 图 像 平 滑 技术 。 噪 声 区 域 
的 大 小 随 着 具体 图 片 而 定 。 通 过 实验 发 现 噪声 区 域 半 径 过 小 会 导致 滤 除 效果 不 好 , 过 大 则 会 
导致 云图 其 它 细 节 受 到 影响 的 情况 。 
1.1.2 时 间 轴 法 : 由 于 云图 数据 是 连续 拍摄 的 视频 帧 ， 每 帧 之 间 间 隔 5 分 钟 左 右 ， 由 此 我 们 
发 现 30 分 钟 内 天 空 灰 度 值 、 星 体 数 量 、 星 光 强 度 等 的 变化 几乎 可 以 忽略 不 计 ， 故 本 文 认 为 
30 分 钟 内 的 云图 所 含 天 空 信息 不 变 。 但 由 于 地 球 自转 和 公转 的 影响 以 及 星星 的 像素 点 较 小 
的 特性 ， 使 得 30 分 钟 内 云图 中 星体 的 位 置 略 有 偏差 。 全 天 相机 的 拍摄 周期 大 约 为 5 分 钟 一 
张 去 图 ， 选 用 与 竺 处理 云图 a 时 间 差 在 30 分 钟 内 的 一 张 云 图 作为 背景 图 b， 裔 历 所 有 像素 
点 的 灰 度 值 ， 发 现 星光 比 背 景 天 空 的 灰 度 值 高 ， 整 体 灰 度 值 越 高 的 图 像 ， 星光 与 天 空 的 灰 度 
值 差 越 大 。 如 果 某 一 点 的 灰 度 值 al, j) 满 足 
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通过 实验 ， 本 文 发 现 Istar 选 取 0.2 去 噪 效果 较 好 。 
(b) 


图 1(a) 星 光 噪声 原 图 (b) 中 值 滤波 结果 图 \c) 时 间 轴 背景 图 (d) 时 间 
Fig.1 (a)Original image of starlight (b) Result image of median filtering (c) Background image of 
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图 b 灰 度 值 平均 值 ，Istar 为 观察 全 部 数据 集 后 
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1.2.1 差分 法 : 差分 法 可 用 于 检测 
期 但 相同 时 刻 
明亮 无 云 
光 噪 声 “。 
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] 光 噪声 同样 本 
目标 。 我 们 将 待 处 
图 ,为 了 避免 混入 
的 图 。 将 RGB 的 原 图 与 背景 图 
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timeline (d) Result image of timeline 
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必要 的 干扰 ， 
并 进行 图 片 相 


1.2.2 边界 阅 值 法 : 月 光 灰 度 值 极 高 ， 但 直接 将 高 灰 度 值 置 为 0， 不 但 会 在 图 片上 留 下 一 个 
突 元 的 黑洞 ， 还 破坏 了 云层 或 天 空 信息 。 而 且 月 光 噪 声 通 常会 伴随 痢 散 射 噪声 出 现 。 散 射 品 
声 是 由 于 月 光 受 天 空中 水 蒸气 等 因素 的 影响 , 发 生 散射 , 在 月 光 附 近 形成 的 面积 较 大 亮度 稍 
暗 的 光 党 。 因 此 仪 将 高 灰 度 值 部 分 置 0 难以 完整 除去 月 光 及 散射 噪声 , 所 以 需要 对 月 光 部 分 
进行 补偿 。 通 过 研究 月 光 边 界 的 云层 信息 ， 发 现 边界 点 的 灰 度 值 与 云图 平均 灰 度 值 成 正比 。 
当 像素 点 灰 度 值 a(ijj) 满 足 
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全 识别 月 光 部 分 ， 过 小 会 将 部 分 天 空 或 云层 误 判 为 
月 光 。 本 文 Imoon 参 考 数 值 为 0.68。 将 灰 度 值 换 成 RGB 三 通道 , 重复 上 述 过 程 


可 得 到 RGB 
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1.2.3 反比 例 线性 变换 法 : 将 RGB 云图 转换 为 灰 度 图 后 ， 我 们 发 现 云层 较 厚 的 部 分 月 光 较 


暗 , 云层 较 薄 的 部 分 月 
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因此 ， 


我 们 可 以 得 出 结论 : 月 光 越 亮 的 部 分 去 噪 后 亮度 变化 越 大 。 因 为 月 光 噪 声 的 灰 度 值 与 云层 厚 
度 成 反比 ， 所 以 去 噪 后 像素 点 的 灰 度 值 与 原始 灰 度 值 成 反比 。 
Pn(ijj) = C — Po(i,j) (1-4) 
式 中 PnGj) 为 去 噪 后 像素 点 的 灰 度 值 ，PoGjj) 为 原始 像素 点 的 灰 度 值 ，C 为 常数 。 为 了 保证 
云图 在 去 除 月 光 噪 声 后 不 突 元 , 月亮 与 天 空 的 连接 处 平滑 ,要 求 连 接 处 的 像素 点 满足 去 噪 前 
后 灰 度 值 相 等 的 条 件 ， 即 
Pn = Po = (am + 255) * Imoon (1-5) 
式 中 Pn 为 去 噪 后 月 亮 与 天 空 连接 处 像素 点 的 灰 度 值 , Po 为 连接 处 像素 点 原始 灰 度 值 ，am 为 
云图 灰 度 值 平 均值 ，Imoon 为 观察 后 选取 的 固定 阔 值 。 
将 式 (1-5) 代 入 (1-4) 中 ， 可 得 出 常数 C 的 取 值 ， 从 而 可 得 处 反比 例 线性 变换 法 的 公式 
Pn(i,j) = 2 * (am + 255) * Imoon —Po(L]) (1-6) 
反比 例 线性 变换 法 先 遍历 所 有 像素 点 , 满足 (1-3) 的 像素 点 确定 为 月 光 噪 声 , 代 入 式 (1-6) 
中 , 其 余 像素 点 不 变 。 


图 2 (a) 差分 法 原 图 (D) 差分 法 背景 图 (c) 差分 法 结果 图 CWA URE (e) 边 界 阔 值 结果 图 (RGB)(f) 时 间 轴 原 图 (g) 时 间 
轴 结 果 图 RGB) (h) 反比 例 线性 变换 原 图 G) 反比 例 线性 变换 结果 图 


Fig.2 (a) Original image of difference method (b) Background image of difference method (c) Result image of difference method (d) 


Original image of boundary threshold (e) Result image of boundary threshold (RGB) (f) Original image of timeline (g) Result image of 
timeline(RGB) (h) Original image of inverse proportional linear transformation (i) Result image of inverse proportional linear 
transformation 
差分 法 的 优点 是 可 以 基本 除去 月 光 噪 声 ， 去 噪 效果 优 于 时 间 轴 法 和 反比 例 线性 变换 法 
缺点 是 会 造成 整体 灰 度 值 降低 ,破坏 云层 信息 ， 且 计算 时 需要 无 云 的 背景 图 ， 不 能 处 理 彩 色 


DH 


图 像 。 因 此 差分 法 适用 于 月 亮 未 被 云层 覆盖 的 情况 。 

边界 阔 值 法 的 优点 是 可 以 在 不 破坏 云层 的 情况 Sea naan 
层 信 息 的 完整 性 ， 还 可 以 处 理 彩色 图 像 ; GEOR AY. Dti FP MEET Hosen 
层 履 盖 的 情况 。 

时 间 轴 法 的 优点 是 可 以 同时 除去 月 光 噪 声 和 星光 噪声 , 可 以 处 理 彩 色 图 像 ; 缺点 是 不 能 
完整 去 除 噪声 ， 需 要 背景 图 。 

反比 例 线性 变换 法 的 优点 是 方法 最 简单 ,无需 背 景 图 ,月亮 区 域 与 天 空 连接 平滑 ,不 会 


使 图 像 突 元 ， 图 像 整体 灰 度 值 不 变 。 反 比例 线性 变化 法 去 噪 效 果 优 于 时 间 轴 法 ， 但 不 如 差分 
法 。 该 方法 不 能 处 理 彩 色 图 像 。 
这 四 种 去 噪 算 法 各 有 利弊 ， 在 实际 应 用 中 ， 要 根据 研究 所 需 、 云 图 中 的 云层 信息 、 云 图 
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数据 库 大 小 等 来 选择 合适 的 云图 
1.3 AXIRA 


RETE. 
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日 光 品 声 的 分 布 、 位 置 、 出 现 方式 、 云 层 对 去 品 算 法 的 影响 等 特点 与 月 光 噪 声 非常 相似 。 


故 可 用 月 光 去 噪 算法 除去 日 光 品 声 。 


光 噪 声 像素 点 更 大 ， 灰 度 值 更 高 ， 很 难 完整 除去 日 


光 噪 声 ， 针 对 这 一 特性 提出 一 种 新 的 去 噪 算 法 
密集 时 间 轴 法 是 在 时 间 轴 法 的 基础 上 改进 的 。 由 于 日 光 像 素 点 较 大 ， 可 选择 30 分 钟 内 


密集 时 间 轴 法 。 


的 5-6 张 云 图 来 实现 日 光 噪 声 的 去 除 。 将 原 图 遍历 像素 点 灰 度 值 与 其 余 背 景 图 进行 比较 ， 当 


满足 下 式 


a(i,j) > (am + 255) * Isun (1-7) 
将 背景 图 像素 点 灰 度 值 赋 给 原 图 。 式 中 Isun 为 观察 全 部 数据 集 后 选取 的 固定 日 光 阔 值 ， 
本 文 Isun 参 考 数值 为 0.68。 也 可 以 将 灰 度 值 换 成 RGB 通道 值 进行 彩 图 去 噪 。 


密集 时 间 轴 法 的 优点 是 能 基本 去 除 噪声 ， 可 以 处 理 彩 色 图 像 。 缺 点 是 需要 多 张 背景 图 ， 


比较 复杂 繁琐 ， 图 像 连接 处 不 光滑 。 该 方法 适合 处 理 大 阳 未 被 云层 覆盖 的 情况 。 


图 3 (a) 密集 时 间 轴 原 图 (b) E 


Fig.3 (a) Original image of dense timeline (b) Result image of dense timeline(RGB) 


1.4 雷电 噪声 


长 时 间 轴 结果 图 (RGB) 


雷电 品 声 通常 出 现在 云层 下 方 ， 因 此 ， 不 存在 被 云层 履 盖 的 情况 。 雷 电 噪 声 呈 树 权 状 ， 


整体 细 长 , 根部 粗 , 枝叶 部 分 细 。 利 用 这 一 特性 , 可 以 采用 形态 学 开 运 算法 来 消除 雷电 噪声 。 
形态 学 开 运 算是 图 像 处 理 的 一 种 常用 方法 , 原理 是 先进 行 形态 学 腐蚀 再 进行 形态 学 膨胀 。 


即 先 消除 图 片 中 的 细小 高 亮 的 目 
像 大 小 不 变 。 


标 物 , 


填补 消除 物体 造成 的 空洞 并 平滑 边界 ,同时 保持 图 


形态 学 开 运 算 可 以 基本 除去 雷 


图 4 (a) 原 图 (b) 结果 图 (RGB) 
Fig.4 (a) 


Original image (b) Result image(RGB) 


处 比较 平滑 ， 是 一 种 比较 理想 的 雷电 去 噪 算法 。 


1.5 反光 噪声 


噪声 ， 保 持 云 层 和 天 空 信息 完整 性 ,雷电 与 天 空 的 连接 


反光 噪声 是 由 于 日 光 或 月 光照 射 在 镜头 上 镜头 反光 造成 的 一 种 噪声 。 所 以 反光 噪声 通常 
伴随 日 光 或 月 光 出 现 。 与 散射 噪声 不 同 , 反光 噪声 不 出 现在 太阳 或 月 亮 周 围 , 其 分 布 不 规律 。 
反光 通常 带 有 颜色 ， 故 可 用 RGB 通道 法 去 噪 ， 此 外 ， 反 光 部 分 的 灰 度 值 比 背 景 天 空 高 ， 还 


可 用 边界 闵 值 法 去 品 。 


1.5. 1RGB 通道 法 : RGB 指 的 是 彩色 图 像 的 红 绿 蓝 三 种 不 同 的 通道 ， 反 光 区 域 的 RGB 三 通 
道 值 大 小 不 同 ， 而 其 余部 分 的 RGB 三 通道 值 相差 无 几 。 将 云图 的 三 通道 分 离 ， 求 出 每 个 通 


道 的 平均 值 。 用 最 高 的 通道 减 去 最 低 的 通道 即 可 得 到 反光 部 分 。 随 后 利用 图 片 减法 , 让 原 图 
减 去 反光 部 分 乘 权 值 系数 ， 可 以 除去 反光 噪声 。 

1.5.2 边界 阐 值 法 : 反光 的 灰 度 值 比 背 景 天 空 高 ， 但 与 部 分 云层 灰 度 值 接近 。 为 保持 云层 信 
息 完 整 ， 需要 对 边界 闵 值 法 进行 改进 。 由 于 反光 噪声 的 像素 点 较 小 ， 在 进行 云图 处 理 时 可 以 
人 为 标记 反光 噪声 部 分 ， 仅 用 边界 阀 值 法 处 理 标 记 部 分 。 


图 5 (a) 原 图 ( 含 标记 ) (b) RGB 通道 法 结果 图 (c) 边界 阔 值 法 结果 图 (RGB) 
Fig.5 (a) Original image(with label) (b) Result image of RGB aisle (c) Result image of boundary threshold (RGB) 

RGB 通道 法 的 优点 是 算法 简单 ， 可 以 基本 除去 反光 噪声 ， 当 权 值 系数 较 大 时 还 可 以 除 
去 部 分 散射 噪声 。 缺 点 时 只 适用 于 整体 灰 度 值 比较 低 的 情况 ， 权 值 系数 很 难 确定 为 最 佳 值 。 

边界 阔 值 法 的 优点 是 能 完整 除去 反光 噪声 ， 不 破坏 天 空 云层 信息 ， 适 用 于 各 种 云图 。 缺 
点 是 非常 繁琐 ， 需 要 进行 标记 。 

1. 6 特殊 噪声 

全 天 相机 云图 中 包含 许多 特殊 噪声 ， 本 文 主要 介绍 其 中 的 四 种 ,包括 积 雪 噪声 、 和 白光 品 
声 、 雨 水 噪声 和 冰晶 噪声 。 

由 于 积 雪 、 白 光 、 雨 水 和 履 盖 整 个 镜头 的 冰晶 会 完全 破坏 云层 和 天 空 信息 ， 故 将 其 作为 
废弃 云图 ， 不 进行 处 理 。 本 文 仅 讨 论 冰晶 噪声 较 小 的 情况 。 

由 于 冰晶 噪声 的 饱和 度 较 高 , 灰 度 值 与 天 空 部 分 相差 无 几 。 可 以 将 饱和 度 和 边界 阔 值 法 
结合 起 来 去 除 冰晶 噪声 。 首 先 定义 饱和 度 阔 值 Su。， 再 结合 边界 阔 值 法 ， 当 像素 点 的 灰 度 值 
KFH Ioe 且 饱和 度 大 于 饱和 度 阔 值 Sige 时 ,判定 其 为 冰晶 噪声 本文 Is 参考 数值 为 150， 
Sice 参考 数值 为 80。 将 冰晶 噪声 的 像素 点 灰 度 值 用 噪声 与 天 空 或 云层 的 边界 点 灰 度 值 替 代 ， 
可 以 去 除 噪声 。 但 云 的 饱和 度 与 冰晶 噪声 类 似 ， 该 方法 会 破坏 云层 信息 ， 故 在 进行 云图 处 理 
a: 时 可 以 标记 冰晶 噪声 部 分 ， 仅 对 标记 处 进行 去 噪 处 理 。 可 以 将 灰 度 值 换 成 RGB 三 通道 ， 重 
e 复 上 述 过 程 ， 得 到 RGB 结果 图 。 
= AURI SEES A FL f TT D c s e fS, 但 连接 处 有 灰 度 值 的 跳 变 , 连接 不 平滑 
© 使 图 像 突 无 。 而 且 该 方法 繁琐 ， 需 要 获取 连接 处 砍 度 值 信息 和 进行 标记 。 


一 


图 6 (a) PR (b) 白光 (WK 〈d) vk CK) (@) 冰 品 ( 小 ) 原 图 (f) 结果 图 (RGB) 
Fig.6 (a) snow (b) white light (c) rain (d)ice (large) (f)Original image of ice (small) (g) Result image (RGB) 
dd 噪声 种 类 及 其 特点 


Table 1 Types of noise and their characteristics 
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噪声 种 类 分 布 区 域 大 小 其 它 
星光 噪声 无 规律 极 小 JEU s 
Hats 靠近 镜头 边缘 ， 位 置 随时 间 BK 伴随 散射 噪声 ， 分 为 有 无 去 
al 层 遮挡 两 种 情况 
光 噪 声 靠近 镜头 边缘 ， 位 置 随时 间 — 比 月 光 噪声 大 伴随 散射 噪声 ， 分 为 有 无 去 
al 层 迹 挡 两 种 情况 
雷电 噪声 云层 下 方 不 定 TABUE, WAN BEHE 
细 
反光 噪声 距 太阳 、 月 亮 有 一 定 距离 较 小 伴随 日 光 、 月 光 出 现 ， 明 彩 
色 
积 雪 噪 声 Bis GE 全 部 不 能 得 到 云层 和 天 空 信息 
白光 噪声 BL REA GE 全 部 不 能 得 到 云层 和 天 空 信息 
雨水 噪声 8 MEA BEL 全 部 不 能 得 到 云层 和 天 空 信息 
冰晶 噪声 不 定 ， 本 文 重点 讨论 冰晶 品 不定 饱和 度 比 背 景 天 空 高 


表 2 不 同 噪声 去 噪 算法 总 结 
Table 2 Summary of different noise denoising algorithms 
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噪声 种 类 去 噪 算法 
星光 噪声 中 值 滤波 法 完整 去 除 星光 噪声 单 便捷 
时 间 轴 法 背景 图 ， 更 好 保留 天 空 云层 信息 完 
月 光 噪 声 差分 法 景 图 ， 适 合 月 亮 未 被 云层 覆盖 情 
WH BELT 亮 被 云层 覆盖 情况 
时 间 轴 法 背景 图 ， 同 时 除去 星光 月 光 噪 声 
反比 例 线性 变换 法 走 ， 连 接 平 滑 
习 光 噪声 差分 法 噪声 
FEY ERS es 
时 间 轴 法 光 噪 声 
反比 例 线性 变换 法 光 噪 声 
密集 时 间 轴 法 需要 多 张 背景 图 ， 基 本 去 除 日 光 噪 声 
雷电 噪声 Base tek 基本 去 除 噪声 
反光 噪声 RGB 通道 法 适用 于 整体 灰 度 值 低 的 云图 ， 简 单 便捷 
边界 闵 值 法 需要 标记 噪声 区 域 
冰晶 噪声 WH BELT 结合 饱和 度 ， 需 要 标记 噪声 区 域 


2 ”交互 式 多 类 型 噪声 去 除 系统 


云图 去 噪 的 算法 较 多 , 需要 根据 实际 情况 选择 不 同 的 去 品 算 法 , 才能 获得 最 佳 去 噪 结果 。 


但 这 样 需要 的 执行 时 间 较 长 , 而 且 部 分 算法 需要 一 张 甚至 多 张 
为 殉 服 以 上 缺点 ， 本 文 将 不 同 噪声 的 去 品 算 法 进行 
量 云图 的 去 噪 工作 。 
通过 对 不 同 噪声 特点 的 研究 , 我 们 发 现 日 光 和 月 光 的 灰 度 值 较 高 , 易于 与 背景 天 空 区 分 


HARS 


繁琐 。 


去 除 系 统 ， 可 以 实现 大 


图 和 噪声 标记 ， 比 较 复 杂 


提出 一 种 交互 式 多 类 型 噪声 


而 其 余 噪声 的 灰 度 值 不 具备 明显 特征 ， 难 以 同 云层 及 天 空 区 分 。 因 此 ， 只 能 判断 云图 中 
否 含有 日 光 、 月 光 ， 很 难 自动 判别 出 星光 、 雷 电 和 反光 噪声 。 
夏 本 文采 用 半自动 的 交互 式 系统 设计 ,通过 人 工交 互 来 判断 是 否 含 某 类 噪声 然后 让 系统 
进行 相应 的 自动 去 噪 处 理 。 
在 这 个 综合 系统 中 , 实际 选择 了 中 值 滤波 法 进行 星光 噪声 去 噪 , 选择 反比 例 线性 变换 法 
进行 月 光 噪声 和 日 光 噪声 去 噪 ， 选 择 形态 学 开 运 算 进行 雷电 噪声 去 噪 ， 选 择 RGB 通道 法 进 


行 反光 噪声 去 品 ， 系 统 流 程 图 如 图 7 Pra. 


T 


是 否 含有 月 光 _ 是 
日 光 噪声 ? 


否 


形态 学 开 运 算 


输出 去 噪 结果 图 


图 7 系统 流程 图 


201909.00185v1 


chinaXiv 


ChinaXiv& (ERAT 


Fig.7 System flow chart 

由 于 部 分 云图 不 含 噪声 ， 并 且 通 过 观察 ， 发 现时 间 间 隔 较 短 的 云图 中 噪声 种 类 、 位 置 、 
大 小 等 信息 相似 ， 为 更 好 的 说 明 该 系统 的 去 品 效 果 ， 本 文 对 原始 数据 进行 了 筛选 。 选 出 200 
张 具有 典型 噪声 的 云 网 , 通过 交互 式 多 类 型 噪声 去 除 系统 的 处 理 , 采用 图 像 质 量 主观 评价 法 ， 
请 20 位 观察 者 对 去 噪 图 像 进行 评价 ， 取 平均 值 作为 图 像 去 噪 效 果 的 结果 ， 如 表 3 所 示 。 其 
中 指标 “去 噪 完整 性 ” 定义 去 噪 后 图 中 无 明显 噪声 的 云图 为 优 ， 去 噪 后 残余 较 少 噪声 的 云 
图 为 展 ， 其 余 为 差 。 指 标 “ 云 层 天 空 信息 完整 性 ”定义 去 噪 后 基本 保持 云层 天 空 信息 完整 
的 云图 为 优 ， 去 噪 后 较 少 云层 或 天 空 信息 被 破坏 的 云图 为 良 ， 其 余 为 差 。 


表 3 交互 式 多 类 型 噪声 去 除 系统 去 噪 效果 
Table 3 Denoising effect of interactive multi-type noise removal system 
去 噪 效果 / 张 去 噪 完整 性 云层 天 空 信息 完整 性 
优 115 141 
良 74 19 
差 11 20 


总 数 200 200 


通过 对 实验 结果 的 “去 噪 完整 性 ”” “云层 天 空 信息 完整 性 ”两 个 方面 进行 分 析 。 我 们 发 
JL. 通过 去 噪 系统 的 云图 多 数 可 以 基本 除去 噪声 ， 去 噪 效果 不 好 的 情况 中 ， 部 分 是 因为 存在 
特殊 噪声 ， 其 余 均 为 反光 噪声 难以 完全 去 除 。 这 说 明 目 前 的 反光 噪声 去 除 方法 还 有 竺 改善 ， 
在 今后 的 研究 中 应 将 如 何 除去 反光 噪声 作为 一 个 重点 。 同 时 ， 该 算法 可 以 基本 保持 云层 和 天 
空 信息 的 完整 性 。 实 验 中 有 20 张 效果 较 差 的 云图 ， 通 过 观察 分 析 发 现 均 为 将 云层 误 识 别 为 
反光 噪声 ， 因 此 如 何 有 效 区 分 反光 噪声 与 云层 仍 需 进 一 步 探索 。 若 要 更 加 有 效 去 除 噪声 ， 应 
该 结合 实际 情况 选择 更 适合 的 去 噪 算法 。 

该 系统 的 缺点 是 不 能 自动 判别 噪声 ， 在 接 下 来 的 工作 中 ， 可 以 考虑 利用 图 像 分 类 技术 ， 
比如 CNN"Y， 进 行 云 图 特征 值 提 取 ， 通 过 距离 分 析 、 聚 类， 实现 噪声 自动 分 类 。 为 提高 分 
类 准确 性 ， 可 用 生成 对 抗 网 络 (GAN) "生成 与 原始 数据 有 相同 分 布 的 样本 ， 增 大 训练 集 ， 
对 系统 进行 预 处 理 。 


本 文 经 过 大 量 的 人 工 研读 ， 将 全 天 相机 云图 的 固有 噪声 进行 了 分 类 ， 共 分 为 星光 噪声 、 
月 光 噪 声 、 日 光 品 声 、 雷 电 品 声 、 反 光 噪 声 和 特殊 噪声 ， 并 分 别 对 其 进行 了 分 析 ， 总 结 出 不 
同 噪声 的 分 布 、 区 域 大 小 和 其 它 特征 ; 然后 根据 这 些 噪 声 的 特性 并 结合 云图 中 云 的 位 置 、 天 
空 信息 等 情况 , 有 针对 性 地 提出 了 8 种 对 应 的 去 噪 算法 , 能 较 好 地 实现 各 种 情况 下 不 同 噪声 
的 去 除 。 同时， 为 了 减少 执行 时 间 ， 快 速 实现 大 量 云图 的 去 噪 ， 本 文 设计 了 一 种 交互 式 多 类 
型 噪声 去 除 系统 , 并 在 实验 中 选用 了 200 张 典型 云图 进行 了 处 理 。 该 系统 的 缺陷 是 不 能 自动 
判别 噪声 类 型 , 在 接 下 来 的 研究 中 需要 考虑 利用 图 像 分 类 技术 实现 噪声 目 动 分 类 ; 而 且 系 统 
除去 反光 噪声 的 部 分 不 是 十 分 理想 ， 需 要 进一步 研究 能 够 有 效 区 分 反光 和 云层 的 方法 。 
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Research on denoising algorithm for camera image 


of all-day camera 


Gao Zhenbin!, Li Mengci’, Qiu Bo!" Cao Zihuang’, Song Tao! 
(1.Hebei University of Technology, Tianjin, 300401 ; 
2. National Astronomical Observatory, Chinese Academy of Sciences, Beijing, 100049) 

Abstract:The all-day camera cloud image has important value in the field of astronomical 
observation, but the noise contained in the cloud image will have an adverse effect. The objective 
of this paper is to perform cloud image denoising processing.This paper has analyzed starlight 
noise, moonlight noise, daylight noise, lightning noise, reflective noise and special noise. And 
summarize the characteristics of different types of noise distribution, region size and so on. Finally, 
eight kinds of image processing denoising algorithms such as boundary threshold method, time 
axis method and inverse proportional linear transformation are proposed based on the relative 
position relationship with the cloud layer. And combine them to form a complete denoising system. 
Through experiments, different algorithms are compared to obtain the best denoising scheme 
under different noise conditions. According to the experimental results, we can draw conclusions 
that the system can effectively remove common noise, basically maintain the information integrity 
of the cloud layer and the sky, and can easily and easily realize the cloud image denoising. 

Key words: All-day camera cloud image; noise analysis; denoising algorithm; image 


processing 


